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Recent studies show that many types of human social activities, from scientific collaborations to 
sexual contacts, can be understood in terms of complex network of interactions. Such networking 
paradigm allows to model many aspects of social behaviour with relatively simple computer models. 
The present work investigates the influence of single leaders on opinion formation within simulations 
of agent based artificial networked societies. Several types of network systems (among them random 
networks, highly clustered, small world and scale-free) are studied. The strength of the social 
influence of individuals is assumed to be dictated by distance from an agent to another, as well as 
individual strengths of the agents. We study the effect of different topologies on the conditions of 
general acceptance of leader's opinion by the society. 



1. MODELING OPINION FORMATION IN 
AGENT BASED COMPUTER MODELS 

Computer simulations are quickly becoming recognized 
and accepted tool in social sciences, viewed as a way of 
testing hypotheses and predictions and revealing some 
simple mechanisms underlying complex behaviour. The 
studies are especially useful when the studied social phe- 
nomena can be quantified, and when enough empirical 
data exists to compare the numerical models with real- 
ity. 

One of such areas is opinion formation in societies, 
measurable both on the large scale (for example through 
election results or polls) or on small scale direct experi- 
ments. Computer simulation models of opinion forming 
date back to 1960s, however, the particular model we 
base our research on is the one of Nowak et al. devel- 
oped further byNowak and Lewenstein Q and Kacperski 
and Holyst Holyst et al. The general approach 

analyses the formation of public opinion through interac- 
tions between individual members of the society, taking 
into account differences in receptiveness, strength of in- 
fluence and preexisting biases. The original work [J has 
shown that interesting macroscopic behaviour can result 
from simple microscopical properties of computer agents. 

The basis for the models in 0, Q, H H Q , which will 
be generally followed here, is the following: 

• A set of N interacting agents form a society. In- 
teractions take form of one to one activities. The 
dynamics of the total system is assumed to take 
place through discrete time steps, with updates 
done globally. 

• Each agent i has, at a given time, his 'opinion' Cj. 
The changes of aggregates of 'opinion' within the 
society or certain subgroups is the topic of the re- 
search. 

• Each agent is characterized by the strength of his 
possible influence on other agents. This allows to 
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model situations of uneven distribution of influence, 
such as presence of strong leader(s). 

• As a special case, one of the agents (the leader) 
is assumed to have the strength of influence much 
greater than the rest of the agents. 

• Interaction between particular agents is governed 
by the strength of agents as well as the remote- 
ness of the interacting agents. We use the term 
'remoteness' in abstract sense, to denote the mea- 
sure of the social separation between members of 
society and its effect on transmission of opinions. 
The overall social impact of agent i on agent j is 
given by combination of acting agent strength and 
the separation between i and j. 

• The agents interact and influence each other in 
turns, changing their opinion after each full turn 
of interactions take place. Agents are allowed 
to interact with themselves, this mimics the phe- 
nomenon of self-support, or inclination to hold 
agent's present opinion. 

• The model allows for extra-social influence or bias, 
applying it uniformly to all agents. 

• The models may allow for the noise in communica- 
tion and changing individual opinion by adding an 
equivalent of temperature to the simulations. 

It is extremely interesting that the results of the works 
cited above have significant dependence on the way the 
'remoteness' between the agents is introduced. The 
works of Nowak et al. Q, Q have used simple spatial 
two-dimensional (2D) model of agent distribution. It has 
been argued that 

'People are more likely to interact with neigh- 
bors, that is those who live close to them in 
physical space' . . . 'Our choice of a 2D lattice 
represents quite well the physical distribution 
of people on flat surfaces. The results of stud- 
ies conducted in Boca Raton, Warsaw and 
Shanghai, have shown that the probability of 
social interactions is decreasing as a square of 
physical distance' ( 2|). 
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Results obtained for the 2D lattice based society have 
shown that there are localized 'bubbles' of uniform opin- 
ion around strong leaders, growing and merging or dimin- 
ishing with changes of external influence. The geometri- 
cal model on a 2D disc has been used by Kacperski and 
Holyst @, with the measure of decrease of the influence 
with distance studied mostly using linear relation. 

In their later works of Kacperski, Holyst and 
Schweitzer 0, Q have argued that the social distance, 
which does not have to fulfill the same conditions as ge- 
ometrical one (such as, e.g. triangle inequality) should 
be modelled by more general and flexible model. The 
authors introduce the notion of social immediacy my 
(between agents i and j ). The immediacy my does 
not need to be equal rriji and the distribution of values 
of my may be arbitrary. In Q three immediacy dis- 
tributions were analysed: uniform random distribution 
> rriij > 2m, exponential distribution and discrete 
multimodal distribution. The results show rapid jumps 
in majority-minority proportion of opinion and presence 
of local clusters of opinion, which in abstract social dis- 
tance space correspond to Nowak's bubbles of opinion. 



2. OPINION FORMATION IN NETWORKED 
SOCIETIES 

2.1. Network systems — an overview 

Recent studies of various networked systems (see for 
example Strogatz Jf}, Albert and Barabasi Q, Dorogovt- 
sev and Mendes U Dorogovtsev et al. [jj, Dorogovtsev 
and Mendes Newman (Til fl^ . Newman and Park 
[lflj). including biological and social systems, the Inter- 
net and World Wide Web have shown several universal 
and interesting effects. The most important among them 
are 'small world effect', degree distributions exhibiting 
power-law or truncated power-law behaviour and signifi- 
cant clustering effects. Large number of results indicate 
that the formation and properties of networked systems 
found in nature — including social networks — shows 
certain regularities and rules. In this paper we study the 
same basic model of opinion formation as Nowak, Holyst, 
Kacperski and co-workers, but we place the agents on 
nodes of several specific networks, differing in topology 
of connections. The different network topologies result in 
significant differences in network characteristics, which 
might influence the process of opinion formation. The 
main networks used in this work are described below. 

In random (RAND) networks the links are distributed 
randomly, thus we have a meshed network of links, with 
agents differing in number of the neighbours, and no gen- 
eral structure. On the average, the maximum distance 
between agents in random network grows as logiV and 
the number of nieghbours per agent is given by the Pois- 
son distribution. 

The nearest neighbor (NN) nets are formed by linking 
together a fixed number of closest neighbors. The easi- 



est way to visualize such network is to place the agents 
on an imaginary circle and connecting each agent to n 
neighbors. Interesting property of NN networks is that 
for small filling factors, agents on the opposite points at 
the circle to communicate must go through many inter- 
mediaries. For 4000 agents and number of neighbors set 
at 10, the longest 'distance' is 200 'hops'. Due to such 
large separation any change of behavior of the agent i is 
seen immediately only by his closest neighbors but the 
majority of the society is affected only after filtering by 
intermediate agents located far from i. 

The small world Strogatz- Watts (SW) networks, intro- 
duced and popularized in recent years, reproduce a curi- 
ous fact observed in many natural and human-produced 
networks, namely that the distance between any two 
nodes of the network, measured in number of links needed 
to connect them, is usually much smaller than that in 
nearest neighbor or even random networks of the same 
filling ratio [14]. The name of the network category comes 
exactly from such observation. One of the ways to real- 
ize a SW network is through a simple reworking of the 
NN model: one cuts a small (even very small!) fraction 
of the links from between nearest neighbors and applies 
them instead between random agents. Keeping in mind 
the visualization of NN networks as nodes along a circle, 
this corresponds to adding connections that criss-cross 
the circle at random. Due to such shortcuts, even if their 
number is very small, the average distance between any 
two nodes drops dramatically. Thus we have a network 
that for each agent, locally is very similar to NN model 
(as most of the neighbors are, in fact, the same), but 
globally the communication through the network is much 
faster. 

The scale-free Albert-Barabasi network (AB) repro- 
duce another effect found in natural and artificial sys- 
tems. In AB networks the dispersion of the number of im- 
mediate neighbors (degree distribution) scales according 
to power law rather than exponential or Gaussian. This 
type of networks results from interplay of two processes: 
growth of the network and preferential attachment (the 
'rich get richer' principle). As a result, the AB networks 
exhibit much greater presence of highly connected agents 
than other types of networks considered here. 

We study also two types of hierarchical networks, sim- 
ple one (HS), with each agent having one link up (to its 
'parent' on higher hierarchy level) and fixed number of 
links down (to descendants), and clustered hierarchical 
network (HCL), in which in addition to links up and 
down, the agents having the same parent are all linked 
together (within the same level). While the simple hi- 
erarchy shows no clustering at all, the clustered version 
preserves the division of the society into separate groups 
and levels, while providing high clustering ratio within 
groups 

For completeness, in addition to the abstract space net- 
works, we study also two traditional spatial geometries, 
2D square array (2DSQ) and 3D cubic array (3DCU). 
In both cases agents may be linked with all immediate 
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neighbors, including diagonal ones (thus the 2D agent can 
have 8 neighbours, and in 3D case up to 26 neighbours). 
We have used periodic boundary conditions. 

In our simulations we have tried to compare results for 
different networks but the same average number of neigh- 
bors Nn- For the regular 2D and 3D networks, where 
the geometrical number of neighbours is fixed, to achieve 
the value of Nn correspondig to other networks, we have 
randomly removed appropriate number of links. The only 
exceptions were the hierarchical structures, where to pre- 
serve the essential top-down asymmetry of the network, 
we have had to accept the lower number of connections at 
the lowest level (edge of the network due to finite size) , 
where only links up (HS) or links up and in horizon- 
tal cluster (HCL) are present. The hierarchical networks 
have therefore smaller average number of links, but pre- 
serve the Nn within the core of the network. 

Figures ^ an d HI compare the two main properties of 
the discussed networks relevant to the present research. 
The results shown were calculated for the average num- 
ber of neighbours Nn = 6. The first figure compares 
the distribution of immediacies calculated according to 
Equation E] (my = 1/dij, i.e. with a = 1) for various 
networks. In all cases one sees very pronounced maxima 
in distribution of my (notice the lin-log scale), but the 
most probable values differ: they are the highest for the 
RAND and AB networks, smaller for SW, 3DCU and HS 
networks, and very small for NN, 2DSQ and HCL. In- 
creasing a would shift the distributions of my to even 
smaller values. It is worth noting that in the NN system 
case due to very large distances most my values are very 
small indeed. 

Examining the distributions of the number of neigh- 
bours one sees that most networks show some limited dis- 
persion around the average Nn = 6, the only exception 
being the scale-free network of Albert-Barabasi. Out of 
the 'normal' networks, the dispersion is largest for RAND 
system, while for NN network the number of neighbors is 
fixed Nn = 6. The AB network contains quite a number 
of agents with high connectivity, even up to 250 in this 
particular case. This turns up as the origin of one of the 
interesting results of our simulations. 



2.2. Use of networks in modelling opinion 
formation 

Even if we concentrate only on studies of strictly social 
phenomena, there is growing body of data that the com- 
munication processes in human societies are best mod- 
elled b y co mplicated network arrangements (for examples 
see |2jll3])- The 'small world' theorem, well established 
experimentally, stating that the separation between ran- 
domly chosen members of society is much smaller than 
expected from geometrical and NN models, and even 
from randomly linked model, can in obvious way influ- 
ence the process of communication between members of 
society, and thus the formation of opinions, or the 'dis- 



tance' from leaders to 'normal' people. On the other 
hand, the presence of extremely highly connected nodes 
(seen in Albert-Barabasi models) gives particular influ- 
ence to some individual members of the society. 

The aim of the current work is to present some results 
of the same type of simulations used before on geomet- 
rical or abstract connections topologies applied here on 
various types of networks. The differences between the 
general network properties might shed some light on the 
way the opinion shift occurs, such as the Nowak's 'bub- 
bles of opinion' in 2D model. 



3. DETAILS OF THE MODEL 

As mentioned above, our model followed almost ex- 
actly that of Holyst and Kacperski. There are N agents 
in the simulated system (we have used A?" = 4096). Each 
agent i is described by strength Si (unchanged during the 
simulation run) and the value of opinion, described by 
<7j = ±1. The process of establishing the opinion in the 
society is modelled in discrete time steps. At each step 
the opinion of any given agent is updated in accordance 
with the combined influence of the other agents, agent's 
own self-influence and possible external conditions. The 
interaction between agents, for example agent j influence 
on agent i is given by combination of agent j strength Sj 
'remoteness' of the agents, described by immediacy my. 
In our work we set a very natural model for my in net- 
worked system 

1 

m u = js (!) 

ay 

where dy is the distance between agents i and j measured 
in number of network links and the exponent a deter- 
mines the ratio of decrease of immediacies with increasing 
distance (we have simulated societies setting a = 1,2,3. 
As ay = 1, 2, 3 . . . the immediacies are always < 1. The 
self-influence terms my were set at 2 for all agents, which 
corresponds to relatively high tendency to preserve one's 
own opinion. 

The strengths Sj were chosen randomly from interval 
[0, 2s]. The only exception was a single agent, called the 
leader, for whom the strength sl was much greater than 
the average s. Moreover, we have fixed <tl = 1 and set 
m-LL ^ 1 so that the leader's opinion <jl would remain 
unchanged during the simulation runs. 

The states of the agents are updated synchronously in 
discrete time steps according to the following formula: 

f 1 with probability cxp(J . / r^ ( !ff_ J . /r) 
a l {t+M) = 1 

[ 1 with probability cxp{I jT)+e£~h/T) 

(2) 

where T is the measure of randomness ('social tempera- 
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ture') and 



N 



h. 



(3) 



The value of h measures uniform external influence on 
the system. 

It is useful to introduce here the rescaled values. The 
scaling is given by 



Nms ' 



T R = 



T 



Nms ' 



(4) 
(5) 
(6) 



with averages m, s excluding the leader. Additionally, we 
define the reference background field B = Nsm, which 
is the maximum value of the background influence of all 
non-leader agents if they all have <jj — 1. Now the social 
impact on agent i, in terms of the rescaled values is given 
by 



h = s R"IhLB + Y sf^-^B + h R B. 
m — ' m 



(7) 



The fact, that through the Equation the network it- 



self establishes values of the immediacies 



allows to 



connect topological properties of the networks (which in 
turn may be closely related to conceptual characteristics 
of the society) to the opinion formation process. Thus, 
our work is much closer to possible experimental and 
practical uses than the previous research. 



4. RESULTS 

In our simulations we were particularly interested in 
searching for conditions under which the leader can con- 
vince a significant part or the whole society to his views. 
We have studied the fraction / of leader's supporters as 
function of unfavorable external influence (h R < 0) and 
leader's strength (rescaled) s R . We have compared re- 
sults of simulations performed with three types of start- 
ing conditions: assuming that initially all agents have 
the same opinion as the leader (oi = 1), random initial 
distribution of Oi and assuming that all agents disagree 
initially with the leader {m = — 1). Results of the three 
types of initial conditions are presented as rows in figures. 

In some of the simulations the leader was chosen ran- 
domly, with the number of immediate neighbors close to 
average, in other cases the highest connected (HC) agent 
was chosen to be the leader. Results are presented in the 
left and right columns of figures EH3 (the only exception 
being the NN network where all agents have exactly the 
same connectivity, so only one set of results is presented) . 
The number of connections of random and HC leader are 
rather similar in most cases — with the exception of the 



scale free AB network, where there may be a difference 
of 1-2 orders of magnitude. 

Figures EH3 present contour plot of fraction of leader 
supporters / for some of the networks studied. The con- 
tours correspond to the following values of the fraction 
of leader's supporters: 0.01, 0.10, 0.25, 0.50, 0.75, 0.90, 
0.99 (from red to green). The social temperature T R was 
set at 0.3. The results, are in most cases remarkably sim- 
ilar. The exceptions, which will be discussed below are 
the NN network and the AB system. 

Figure Ul presents / as function of s R for a fixed value 
of h R — —1.5, for a set of network types and initial 
conditions. Left column corresponds to extremely dis- 



favourable initial conditions (Jo 



T, right column to 



a random distribution with average a = 0. Two sets of 
lines correspond to leaders placed at random (red) and 
in a HC position (green). 

Let's discuss first a typical situation, such as the one 
found in RAND and SW networks, as well as the AB 
network when the leader is chosen randomly. For low 
s R and significant negative h there are no supporters 
of leader's opinion. Increasing s R leads at first to estab- 
lishment of small cluster of supporters (manifested in the 
area where 0.01 < / < 0.25. In figure |H1 this is clearly 
visible as plateaus of /. Further increase leads to rapid 
transition of the whole system to polarized, supporting 
state of / w 1. The transition region is approximately 
given by equation 



s% = -h R + C, 



(«) 



where the value of the constant C depends on the tem- 
perature and initial conditions. This equation simply re- 
flects the strength of the leader needed to overcome the 
average background field of initial configuration. In low 
T R limit C ~ 1 for the <7j = — 1 starting condition and 
C = — 1 for the opposite case; for random starting con- 
dition C — 0. These results are quite straightforward: 
when the leader's strength overcomes the combined effect 
of the background influence and external conditions, the 
agents start to change opinion. With the increase of their 
number the conditions become even more favourable and 
transition to supportive state occurs. 

That the range of values for which the local cluster of 
supporters is present is greater for the favourable starting 
conditions than in the random or disfavourable case is 
also quite evident: due to relatively large value of self- 
influence terms ma in Equation Q it is more difficult 
for the leader to convert an agent than to maintain the 
supporting agent's attitude. 

The results are presented for moderate values of h R 
and s R , comparable to 1. For much larger values the ex- 
ternal influence and the leader's influence dominate any 
local statistical fluctuations due to quenched disorder. 
The dynamics of the system is then trivial, given, for any 
agent i by competition between the leader term s^^t^t 
and h R . Due to discrete nature of distances d iL and 
therefore miL, agents with the same are convinced 
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(for given h R ) at discrete values of s R , and / shows sim- 
ple multi-step behaviour. 

We return now to the special cases. We'll start with 
the NN network. Due to the presence of very large sepa- 
rations between agents in NN system, leader's influence 
on distant agents is extremely weak, and the only prac- 
tical way to supportive state is the gradual 'conviction' 
of the growing range of neighbours. As a result we see 
that the region of / grows linearly with increasing s R . 
The growth of / is slow, and very large values of s R are 
needed to ensure the wholly supportive state. The NN 
society may be called strongly localized, and the presence 
of very high separation between agents clearly influences 
the opinion formation. It should be noted that our NN 
network is, in fact, a 1-dimensional (ID) network. The 
2D and 3D networks, for which the average separation 
grows not linearly with N but rather as ViV and \/~N 
respectively should weaken the effect, but still preserve 
it. Our simulation size of 4096 agents is, however, too 
small to make it visible, as the furthest distance on 2D 
network is 32 and for 3D case only 8. On the other hand, 
for random networks, where the largest distance grows 
as log TV and for SW and AB networks where it is even 
smaller, the conditions for rapid transition are present 
and the such transition is observed. 

The second interesting case is the AB network with the 
leader in high connectivity position. Such leader has a 
lot of neighbours with dy = 1 and thus with my = 1. 
The influence on those agents is relatively strong, and as 
these agents form a significant fraction of the society's 
size (in our case 250/4096 ~ 6%) it is much easier to 
achieve the supportive state 

s* L >K(-h R + C), (9) 

with K < 1. 

Figure [3 presents results for the AB network with 



a = 2. The decrease of my with distance being much 
faster, there are significant differences to the a = 1 case. 
The transition region is here more gradual, with visible 
spread of the contours especially for large negative h R 
region, also for the leader in random position. While the 
contour / = 0.90 falls into similar position as in the a = 1 
case, the / = 0.99 contour is shifted to much higher val- 
ues of s R in the case of disfavourable and random initial 
conditions. This is due to the difficulty of converting even 
moderately remote agents due to rapid decrease of my . 

The results presented in this work indicate that simple 
models of networked societies, based on network types ac- 
tually found in nature and created by humans, together 
with a simple formula for turning social distance to influ- 
ence strength one is able to derive the picture of the con- 
ditions necessary to bring the consensus of opinion due 
to the influence of a single strong individual. For most 
of the networks the simulated behaviour is quite sim- 
ple, but for some of the network topologies (especially 
for the scale-free Albert-Barabasi networks) the results 
are significantly different. The conditions necessary to 
convert the population to leader's opinion depend cru- 
cially on leader's position and connectivity in the net- 
work. Also, for strongly localized nearest neighbour net- 
works the achievement of fully supportive state requires 
far greater leader strengths than for small-world or even 
random networks. 

As the Strogatz- Watts and Albert-Barabasi networks 
are found to be relevant in many human societies and 
activities (such as friendship and acquaintance networks, 
sexual contacts, scientific collaborations, scientific cita- 
tions, internet WEB pages links and even telephone calls 
structure) the presented model may be found of impor- 
tance in the analysis of the processes of opinion formation 
in some of these cases. 
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Figure 1: Distribution of rriij for various types of networks 
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Figure 2: Distribution of the number of neighbours for various types of networks 
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Figure 3: NN network. Fraction of leader supporters / compared to society size, as function of the rescaled leader strength s L 
and external conditions h R . Contours correspond to / = 0.01, 0.10, 0.25, 0.50, 0.75, 0.90, 0.99. T = 0.3 and a = 1. 
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Figure 4: RAND network. Fraction of leader supporters / compared to society size, as function of the rescaled leader strength 
sf and external conditions h R . Contours correspond to / = 0.01, 0.10, 0.25, 0.50, 0.75, 0.90, 0.99. T = 0.3 and a = 1. 
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Figure 7: AB network. Fraction of leader supporters / compared to society size, as function of the rescaled leader strength Sl 
and external conditions h R . Contours correspond to / = 0.01, 0.10, 0.25, 0.50, 0.75, 0.90, 0.99. T = 0.3 and a = 2. 
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Figure 8: AB network. Fraction of leader supporters / compared to society size, as function of the rescaled leader strength s L 
for fixed external conditions h R = —1.5. Left column corresponds to <jj = —1 starting conditions, right column to random <jj, 
with a = 0. Red lines are the results of simulations with leader in a random position, while green lines correspond to leader as 
highest connectivity (HC) position. 



